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Abstract 

We consider electron transfer from a redox to an electrode through and adsorbed 
intermediate. The formalism is developed to cover all regimes of coverage factor, from 
lone adsorbate to monolayer regime. The randomness in the distribution of adsorbates 
is handled using coherent potential approximation. We give current-overpotential pro- 
file for all coverage regimes. We explictly analyse the low and high coverage regimes by 
supplementing with DOS profile for adsorbate in both weakly coupled and strongly cou- 
pled sector. The prominence of bonding and anti-bonding states in the strongly coupled 
adsorbates at low coverage gives rise to saddle point behaviour in current-overpotential 
profile. We were able to recover the marcus inverted region at low coverage and the 
traditional direct electron transfer behaviour at high coverage. 



1 Introduction 



A proper understanding of electron transfer reaction through an adsorbate intermediate 
constitutes the first step towards modelling the charge transfer across a chemically modified 
electrode [1-3], through a molecular wire [4,5], or the phenomenon of the molecular elec- 
tronics [6-10]. In fact the indirect heterogeneous electron transfer is a recurring feature in 
all these processes. 

In the present communication, we consider the kinetics of an adsorbate mediated electron 
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transfer reaction. The adsorbate is taken to be a metal ion. The reactant is supposed to 
couple with the adsorbate orbital alone; the direct couphng between the reactant and Bloch 
states in the metal electrode is neglected. In the present study, the adsorbate coverage 
factor 9 is allowed to take any arbitrary value in the range (0, 1). Thus starting from a single 
adsorbate case, corresponding to 9 — > limit, the formalism remains valid all the way up to 
a monolayer regime = 1) . For metallic adsorbates, the adsorbate orbitals remain spatially 
localized in the low coverage regime. But in the monolayer regime, one obtains extended 
electron states in the adlayer. These states now form a two-dimensional band [11,12]. The 
localized adsorbate states interact strongly with the solvent polarization modes. On the 
other hand, the interaction of extended electron states with polarization modes are much 
weaker, and as a first approximation, it can be neglected [13] . This progressive desolvation 
of adspecies, when the coverage is varied from zero to one, changes adsorbate orbital energy 
by a few electron volts and hence must leave very significant effects on the electron kinetics. 
In addition, in the monolayer regime, the metallic adlayer acts as the electrode surface. 
As a consequence, the adsorbate mediated electron transfer exhibits the characteristics of a 
direct heterogeneous reaction. We investigate how the desolvation and metallization of the 
adsorbate layer influences the charge transfer kinetics. 

The crucial difference between the heterogeneous electron transfer reaction through an 
adsorbate and a direct electron transfer to an electrode arises due to (i) a possible change 
in the electronic coupling term between the participating orbitals, and (ii) modification 
in the relevant density of state. The electronic coupling strength can change due to the 
particular symmetry configuration of the orbitals. Besides, the equilibrium distance between 
the reactants may vary in both the situations [ ]. Next, the density of states (dos) of a 
metaUic electrode is usually broad, and a slowly varying function of energy [ ] . This feature 
enables one to replace the energy dependent dos by its value at the Fermi level. As the 
coupling between the band states in the electrode and solvent polarization modes are usually 
neglected, the electrode dos does not exhibit any temperature dependence. The adsorbates, 
on the other hand have a narrower dos, which depends on temperature. This follows from 
the solvent induced broadening of the adsorbate orbital. When this broadening mechanism is 
absent, adsorbate dos ceases to be temperature dependent. In addition to the solvent induced 
broadening, the adsorbate dos acquires an additional temperature independent width due 
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to the hybridization of its orbital with the Bloch states in the electrode. The location of 
adsorbate dos vis a vis the dos of redox couple play an important role in determining the 
charge transfer kinetics. In fact catalytic effect can be observed when strong overlap occurs 
between these two density of states. 

The adsorbates exhibit different structural arrangements at different coverage. Even 

at a fixed coverage, more than one kind of distribution pattern can be observed in the 
adlaycr [14-16]. Modelling each configuration separately poses a difficult task. Therefore we 
consider a random distribution of the adsorbates in a two dimensional layer. Subsequently, an 
'effective-medium' description is used for the adlayer. This procedure captures the essential 
features of the adlayer in an average sense. 

The plan of the paper is as follows: In section 2, we present the model Hamiltonian and 
the expression for anodic and cathodic current in terms of system parameters, whose detailed 
calculations as shown in appendix. In section 3, the results of numerical analysis is presented along 
with the various DOS for different regimes and also the profile of current at different coverages 
are considered along with explanation for the observed behaviour. In section 4, we summarize our 
results and an overview of the whole work is given 

2 System Hamiltonian and Current 

An adsorbate has strong electronic coupling with the substrate band states as well as 
it has electronic overlap with the neighbouring adspccics. The latter coupling leads to a 
two-dimensional band formation in the adlayer at higher coverage. The solvent-adsorbate 
interaction and surface plasmon-adsorbate interaction, both modelled within the harmonic 
approximation, are responsible for the solvation and image energy for the adsorbate, respec- 
tively. Similar interactions are present for the redox-couple, which is supposed to interact 
weakly with the adsorbate orbital. Taking into account various system components and in- 
teractions among them, an effective Hamiltonian for an adsorbate mediated electron transfer 
reaction can be written as 



^ = Y1 ^^^^^ + + 





a 



a 
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The redox species is coupled to an adsorbate located at a site i — ain the adlayer. Hg is the 
Hamiltonian for the 'electrode - adsorbate - solvent' subsystem 

(2) 

{i} specify sites in the adlayer. k and r label electrode and reactant electronic states, a is 
the spin index and e is the energy value, n, and c respectively denote number, creation and 
annihilation operators for electrons, = 1, 2, 3, 4 label oscillator modes corresponding to the 
orientational, vibrational, electronic solvent polarization and surface plasmons, respectively, 
and uji, are the associated frequencies. 6, 6^ are the annihilation and creation operators for 
the boson modes, v represent the coupling term between the electronic states. A signify 
the strength of adsorbate and reactant coupling with the boson modes. Subscripts o and c 
respectively refer to the reactant and adsorbate core. 

er{bu + bl} = e° + Y^ \r.{b, + bl) (3) 
iia = eaa + ^ Kv {bv + bl) (4) 

V 

e° and are the reactant and adsorbate orbital energies in the gas phase The expression (4) 
gives the energy of adsorbate at site i when it is occupied. In case no adsorbate occupies the 
site i, the expectation value 

< > — > oo (5) 

which ensures no charge transfer through an unoccupied site. While evaluating the shift in 
adsorbate orbital due to its coupling to boson, the boson mediated interaction between dif- 
ferent sites are neglected. Next, since the adsorption of a single type of species is considered, 
we replace Ajj, by \av and Xid, by Xcu- The randomness associated with the site energy can 
be handled using the coherent potential approximation [ ] . 

Treating the magnitude of Var to be a small quantity, the anodic current contribution is 
obtained within the linear response formalism as {cf .Appendix A) 
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/oo 
sgniX^ie, 9)) (1 - /(e)) pTie) p^{e)de (6) 
-oo 

where /(e) = {l + exp{—(3e))~^ is the Fermi distribution function. Here zero of the energy 
scale is set to be at e/ for a direct electrochemical electron transfer reaction. p^{e) and and 
p^"(e) are the adsorbate and the reactant density of states. 

PT{e) - -^Re{w{z)) (7) 



w{z) — e ^ erfc{—iz) (8) 



Z = (-g-" + i|X2(e,e)|)/(2VP), (10) 

_ (.-.o + 5]^)i?;:^ 

= X,(e, e)-el + Y,^ ^—^ (11) 

Re{w{z)) denotes the real part of [w(z)]. Xi and X2 is obtained through the relation (cf. 
Appendix) 

X^{e,e) +1 sgn{X2{e,e))\X2{e,e)\ = Gj^ + K,{e,< >] (12) 
with being the coherent potential which has to be estimated self-consistently. 

— ^cv + ^ov + ^rv (13) 

in case of anodic current. 

V V V 

are the reorganization energy for the reactant, adsorbate, and the cross reorganization energy, 
respectively. 
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Alternatively, we can also write 



el^e'^-J2^^F^^-FS-E;^rj-E; (16) 

V 

where 

F^ = e«-^^-2^^^ (17) 

Fq and Fr denote the free energies of the redox-couple in the oxidized and reduced states. 

(R — (Q = f^r' "^Ri^ — ^ru — ^ou- Thus Fr ~ Fq givcs the overpotential r/ for the electron transfer 
reaction. Similarly, the fraction of overpotential drop between the electrode and adsorbate is related 
to the change in the adsorbate free energy during the reaction 

^^ = e'^^-Y^^=Fl,-F-0-K^ocr,-K^Kr (19) 



Rewriting the Anodic current expression in terms of overpotential, the expression for Q and pr 
takes the form as shown below. 



2e: 



AttEj;: 



exp 



AE; 



(20) 



(21) 



Proceeding along similar lines of argument for the cathodic current, and noting that \i, 
Aci/ + Aoi^ + Xav for cathodic current, the expression for Q and pr obtained as shown below, 



I. 



/oo 
sgn{X2ie, 0)) /(e) p„(e) p,(e)de (22) 
■oo 

e-v-E; + EDEl, 
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2e: 



(23) 



The coupling constants between adsorbate and various oscillator modes are scaled by a factor 
^(1 — 9^) to take into account the disolvation effect as adlayer itself exhibits metallic properties in 
the higher coverage regime. Consequently, the solvation and reorganization energy for the adsorbate 
get scaled by a factor (1 — 0^), and the solvent induced cross energy terms are scaled as 1^(1 — 0'^) []. 
No such scaling is present for solvation and reorganization energies of the redox-couple. Thus the 
scaling laws for the various re-organisation are as follows 

KriO) = VO^KM ; Kio) = {I - 9')e:{o) (25) 
3 Numerical Results and Discussions 



The basic concern in the article is toward currcnt-ovcrpotential characteristics with specific 
emphasis on the variation with the coverage factor (9) and the fraction of overpotential drop {arj) 
across the adsorbate. A first look at the expression for anodic current for a shows that the current 
is an overlap integral of three terms corresponding to the availability of vacant energy level at 
the electrode (1 — /(e)), the density of states of the solvated redox couple p^" and the density 
of states of the adsorbate The redox density of states has a Gaussian form in terms of e. 

The self-consistent evaluation of the coherent potential ka-{9) enforces a numerical derivation of the 
adsorbate density of states. However in the following limiting cases, kcr{9) takes the value 

lim ka = e Wii (26) 

and 

lim ka = Caa (27) 
6>-»l 

where 

Wn = Y.^ (28) 

k ^ 

Consequently, the adsorbate density of states can be analytically obtained in the limits 9 
and 1. Additionally, e^^ involved in performing the self-consistent evaluation of the coherent 
potential takes the value as ar] — Ea{9) + Ear{9) for anodic current evaluation and ar] + Ea{9) 
for cathodic current estimation. 



-P 



(24) 
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In what follows, we describe the current vs overpotential profile for different sets of parameters. 
The adsorbate-electrode interaction is treated both in the weak {v = O.beV) and strong (v = 2.0eV) 
coupling limits. When the coverage is low, the adsorbate density of states has a single peak Fig. [1] 
. An important consequence of the strong coupling limit is the splitting of the adsorbate level in 
bonding and anti-bonding states for low 9 Figl2j This feature is recaptured in the present analysis 
since energy dependence of A(e) is explicitly treated in the present approach (Eqs IA-141 IA-15p . On 
the other hand, the well known wide-band approximation for A(e) fails to provide the bonding anti- 
bonding splitting. In the monolayer regime, due to the 2-d bond formation by the adsorbate layer, 
its density of states acquires a flat profile, irrespective of the strength of the electrode-adsorbate 
coupling (Fig [H [2]) . The table I summarizes the values of parameters used in the calculations. 



Table 1: Values of parameters used in calculation in eV 





V 


A|| 








Ear (0) 


Ea (0) 


strong 


2.0 


0.75 


1.5 


4.5 


1.0 


0.25 


0.75 


weak 


0.5 


0.75 


1.5 


4.5 


0.6 


0.2 


0.4 



Ideally, under zero overpotential condition, the anodic and cathodic currents are expected to be 
equal in magnitude. This implies that the profile of the product /o""'(e) * Pa^i^) for anodic current 
is identical to the product profile Pr^^{e) * /Oa"*(e) for the cathodic current. This is a consequence of 
equal separation between the peak positions of adsorbate and reactant density of states for anodic 
and cathodic processes during equlibrium. [Fig. Oil]. The corresponding plots for strongly coupled 
regime is also shown in Fig. [5l [6] 

As noted earlier, the electrochemical potential p has been set as the zero of energy scale for 
the direct electron transfer reaction. The presence of additional charge particles for the bridge 
assisted electron transfer reaction, namely the adsorbates, changes the equilibrium potential of 
the electrode. This is turn gets reflected as a ^ dependent variation A(j){9) in fj.{= 0). The fact 
that the anodic and cathodic currents at equilibrium potential are identical in magnitude provides 
a novel method for the determination of A(/)(0). Thus the relation Ia{r] = 0) = Idr] = 0) with 
/(e) = {1 + eyLp{-P{e + A(j){9)))-^ { } (cf eqs. Eland [22|) enables us to evaluate A(t>{9). The 
variation of A(j) with respect to 9 is shown in FigjT] in the limit of weak and strong adsorbate- 
electrode interaction, with E; = 0.6 eV, E^^^O) = 0.2 eV, E^ (0) = 0.4 eV. The value of A(j){9) 
depends on the strength of coupling v; its magnitude increases as the coupling becomes stronger. 
|A(/>(0)| is again large for low 9 values and remains almost constant in this region. Note that in 
this regime, the charge on the adsorbate remains localized on the adsorption site. |A(/>(0)| starts 
diminishing sharply for 9 > 0.6 and it tends to as ^ 1. This behaviour is expected. As 9 ^ 1, 
the adsorbate layer becomes metallic and gets incorporated in the electrode. The electron transfer 
acquires the characteristics of a direct heterogeneous reaction, and consequently as noted earlier, 
the electrochemical potential fi again lies at the zero of the energy scale. 

We first present the current-overpotential profile in the weak coupling limit {v = 0.5 eV) for a 
range of 9 and a. The employed values of various reorganization energies are E^ = 0.6, -^^^(O) = 
0.2, E^ (0) = 0.4. The general behaviour can be analysed by looking at the case of lower coverage 
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and high coverage regimes respectively, and then by investigating the effect of variation of a in these 
hmits. Fig. [8] shows that for a fixed a, anodic current as weh as the current peak height increases 
with 9 in the small 9 range (curve a and b). This feature arises due to a better overlap between 
the reactant and adsorbate density of states, whose peak positions are approximately separated 
by a distance E^. + E'^{9) — E'^j.{6). An increase in 6 reduces E^ and E^^ (cf eq. [25]) . and hence 
the peak separation diminishes and the overlap gets enhanced. The presence of anodic current 
peak at r]p signifies negative differential resistance for r] > rjp. This feature is absent in the higher 
coverage limit. For large value of 9, the current at higher r] exhibits a saturation effect. This is a 
consequence of the fact that the maximum n in the adsorbate density of states p^" is now absent. 
yoJJ" now acquires a plateau profile (Fig[T]). The plateau height, and therefore the overlap between 
the reactant and adsorbate density of states decreases with the increasing coverage . Therefore a 
decrease in the saturation current results as — > 1 (curve 9 = 0.7 and 0.9 in Fig. [8|). 

The effect of the a variation on the anodic current is highlighted in Fig. [9l [101 and [UlThe effect 
is more pronounced in the low coverage regime due to the presence of adsorbate density of states 
peak. The reactant and adsorbate density of states peak separation increases with the increasing 
a. Consequently, the maximum overlap between the two occurs at larger r]. This explains the 
occurrence of the anodic current peak at higher r] values as a increases . On the other hand, the 
near constant adsorbate density of states for large 9 ensures a minimal effect of a variation on the 
anodic current (Fig. [TTl [T2]) . 

Next the strong coupling limit with [v = 2.0 eV, E; = 1.0, E^^^.{0) = 0.25, E^^ (0) = 0.75 ]is 
considered. Figures [13|14|15 l and 1161 shows the current overpotential response in the strong coupling 
regime. As in the case of low coverage, the la vrs r] plot exhibits a negative-differentail region (Fig. 

da El). 

More importantly, the presence of two peaks in /o"" when coupling v is large and 9 is small (Fig. 
121) leads to a saddle point and a maximum in the la vrs. ry plot. For the set of parameters currently 
employed, the lA,Max now occurs at a much larger ry in comparison to the weak coupling limit, 
and may not be accessible experimentally. However, the saddle point in the current appears in an 
overpotential range where the anodic current peak appears in the weak coupling limit. For large 
coverage, current potential profile are similar in strong and weak coupling limit. Interestingly, the 
saturation current is smaller in the large coupling case due to a decrease in the height of p^"". In 
fact this lowering of the current in the strong coupling is holds true for any coverage and rj. This is 
shown in Fig. [T71 wherein the variation of equilibrium current 1° with respect to coverage is plotted 
The 1° is smaller for larger v, and as explained earlier in the context of Fig. [8l shows a maximum 
in the intermediate coverage regime. However it may be noted that when — > 0, current would be 
proportional to and an increase in v in this very weak coupling limit will lead to an increase 
in the current. 

The results presented till now correspond to anodic current. But the formalism developed here 
also yields the cathodic current. In fact the equivalence of anodic and cathodic currents at the 
equilibrium potential has been earlier employed to determine the the variation in the equilibrium 
potential due to varying adsorbate coverage. The dependence of the cathodic current Ic on over- 
potential 7] is plotted in Fig. p^ . It is often presumed that Ic = e~'^Ia (Fig. [T8l) . The present 
'microscopic' calculations show that it is not entirely true. The calculated current value is slightly 
larger than the e~'^Ia when 77 < 0. 
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The high coverage regime oi 9 ^ 1 corresponding to a formation of monolayer of A decrease 
in the current for higher rj when the coverage is low virtually mimics the Marcus inverted region 
for a homogeneous electron transfer reaction. On the other hand, the current getting saturated at 
higher r] when the coverage is large is also true for a direct heterogeneous electron transfer reaction. 
Thus depending on the extent of coverage, an adsorbate mediated electron transfer at an electrode 
exhibits the characteristics of both homogeneous and heterogeneous electron transfer reactions. 
The localization of adsorbate electron at low coverage and its delocalization at high coverage is the 
reason behind this phenomena. 

4 Summary and Conclusions 



In this work, we considered electron transfer in an electrochemical system, from a solvated 
redox to an electrode mediate by intervening adsorbate atoms. Further randomness is introduced 
in the model in terms of the coverage factor which relates to the number of adsorbate atoms 
adsorbed on the electrode surface. The theory developed is valid for a range of regime, lone 
adsorbate mediate transfer to the monolayer formated direct electron transfer regime. The inherent 
randomness involved in the adsorbate distribution on the surface has been tackled by coherent 
potential approximation (CPA) and separate expression are derived for anodic and cathodic current. 

Explicit attention was paid to the low coverage and high coverage regime, even though the 
formalism is valid for all regime, since at these two regions the theory could be compared with pre- 
existing literature. Plots were also provided for intermediate regimes and additionally, the effect 
of the adsorbed atoms on the Fermi level of the electrode were incorporated by means of a shifted 
potential A(f){9), ensuring that the anodic and cathodic current were equal under zero overpotential 
condition. 



The analysis also provides a novel method for determining the variation in A(j){9) with changing 
adsorbate coverage. 

The fraction of overpotential drop across the electrode-adsorbate is incorporated and the collec- 
tive plots are analysed. We have proved that this fraction of overpotential drop plays a significant 
role in determining the response behaviour of current, typically the location and extent of the 
maximas in case of lower coverage situations, while in case of high coverage regime, the effect is 
not profound and the electron transfer follows the traditional direct electron transfer as expected 
from heuristic arguments. 

The dependence of anodic current in the weak and strong electrode-adsorbate coupling is ana- 
lyzed. In the former case, vrs overpotential profile exhibits a peak, where as in the later case, 
and in the same overpotential region, the current plot shows a saddle point behaviour. This fact 
can be used to distinguish a weakly chemisorbed bridge from a strongly chemisorbed one. These 
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distinguishing features occur only when the coverage is low. At high coverage, la V plots have 
identical profile for weak and strong coupling cases 

The calculated cathodic current gives a slightly higher value of Ic in comparison to a presumed 
Ic which equals e~'^Ia- 

At low coverage, it is possible to recover the Marcus inverted region, which is absent when the 
coverage is large. The localized nature of the adsorbate orbital when coverage is low, and its getting 
delocalised for high coverages leads to this behaviour. 
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appendix 



The microscopic current associated with the electron transfer reaction depends on the average 
value of the rate of change of electronic occupancy of the redox orbital [ ] 

Treating small quantity, a linear response formalism can be used to evaluate < hr >. 

Consequently, 

I = IflY. < [^/UO)' ^^-(*)]- > (A-2) 

n ^ J-oo 

where 

Via = VarClCr (A-3) 

The first term in the commutator leads to anodic and the second one gives the cathodic current. 
The expectation value in (A-2) now corresponds to a density matrix defined for the Hamiltonian 
H' = H — ^^{V^^. + Vjr)- H' also determines the time evolution of various operators in (A-2). 
Employing the Frank-Condon approximation and treating the low frequency polarization modes in 
the semi-classical approximation, anodic current is obtained as 

la = T dt\Vat\'' « Cl{0)Crit) >f< Caa{0)Cl{t) >F>B (A-4) 

" J -oo 

where 

1 poo 

< cl{0)criT) >F= - e-''^'H{e - er)de (A-5) 

TT J-oo 

The time correlation function involving c^, cl can be expressed in terms of adsorbate Green's 
function 

< Caa{0)cUt) >F= - (1 - f{e))e''^l\lmGii)i=ade (A-6) 

^ J-oo 

where 

{Gii{e))i=a = < ^\Caa < 777 >c,i=a cL|0 >F (A-7) 

e — rl 

Here < ... >f implies an average over electronic degrees of freedom, keeping bosonic variables as 
fixed parameters. < ... >b denotes the thermal average over boson modes. < ... >c,i=a denotes 
a restricted configuration average. It implies that while obtaining the configuration average, the 
site a, which is occupied by an adsorbate and through which the electron transfer takes place, is 
excluded from the averaging. The occupancy status of the remaining sites are still unspecified. We 
replace the random medium encompassing the remaining sites by an effective medium using the 
CPA technique. The picture which now emerges is the one in which a reactant is coupled to an 
adsorbate occupying the site a, and this particular adsorbate is embedded in a two dimensional 
effective medium. 



12 



The randomness associated with these sites can be handled using the coherent potential ap- 
proximation [ ]. Accordingly, the random energy operator eia{< Qu >)nna in eq. |4]is now replaced 
by a deterministic operator k^rii^. The coherent potential /^^(e, < >) is same for all the sites, 
but depends on the energy variable e. is determined self-consistently through the expression [ ] 

Q.. ^ J_Y \ = (A-8) 

iV|| ^ e - k„{e, <qy>)-eu- W'{e, u) e^a - k„{e, < >) 

where 2D adsorbate lattice has A^|| number of sites, and 

W{e,u) = EJe'''■^Hv^J+W|^{e)] 

= eu + Eje'^-^'-Wl^i^) (A-9) 

= e^ + W'{e,u) 

The {Gii)i=a can be related to the complete configuration averaged GF Gu 

{G^i)^=a = Gu{e)[l - {e^M - kA<^))Gii{e)]'' (A-10) 

The tedious summation over the momentum k of metal states and the momentum u of the 
Bloch states in 2D adsorbate layer, commensurate with the underlying electrode surface lattice, 
can be considerably simplified under the following assumptions, (i) The separability of the metal 
state energy in the direction parallel and perpendicular to the surface, (ii) The substrate density 
of states in the direction perpendicular to surface is taken to be Lorentzian, whereas the same is 
assumed to be rectangular along the surface, (iii) The adsorbate occupies the 'on-top' position 
on the electrode, and is predominantly coupled to the underlying substrate atom with coupling 
strength v. Consequently, eq.(A9) now becomes 



Gi 



1 



eaa{< Qu >} - ka{e, < qu >} 



1 (^-C)ln(^)-(5-C)ln(|^ 



2A||(5-^);U 



(A-11) 



with 



A/B = h{C + D)±{{G + Df -A{GD--)Y2] (A-12) 



C = ' ^ ; D = e-iAr, M = 96/ Au (A-13) 

5 is the half-bandwidth of adsorbate monolayer and fi is the same for an arbitrary coverage 
e. 2A|| is the sub 
2(A|| + Aj_) = At. 



2A|| is the substrate bandwidth at the surface and A^ the total bandwidth of the substrate 



1/2A|| ; -e < All < e (^"1^) 
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Where the separabihty of metal states energies in directions parallel and perpendicular to 
the surface is a crucial assumption made here, = + e^^ . 

The self-consistent value of k(j(e) is obtained from eg. I A- 111 which also determines Gu. 
Expressing 

GT^^ + K^{e,<q,>} = Xi{e,e) +i sgn{X2ie,e))\X2{e,9)\ (A-16) 
enables us to write (cf. eq.(All)) 

Im{Gu)^=a = -sgn{X2) ^^'\ (A-17) 

[Xi{€,e) - eaa{qu)) +X^ 

where sgn(X) = +1 when X > and -1 otherwise. 

To evaluate the anodic current /q, one finally needs to carry out the thermal average over low 
frequency boson modes. The required density matrix for this average in the semi-classical limit is 

/oo 

W{q,)dq^ (A-18) 
-co 

W{q,) = exp[-l3 '^{Pi + <11)+Kqu\ (A-19) 

!^=l,2 



With the above defined probability function the net expression for anodic current is shown 
below 

la = e{nr) I Var P j de j (U^dq^) j dtj dT[{l - /(e)] 

exp[-P ^ Y + + ^"'i"^ exp[i{e - - ^ Kuqu)t] 

u=l,2 V 

exp\i(X^{e,e)-ea-Y.\avqv)r-\X2{e,Q) || r |] (A-20) 



where Z = J {JV^dq^^ex-py— 15 ^{Pu + Qu) + ^i^Qiy] The expression for cathodic current has a 
similar form with the 1 — /(e) replaced by /(e) and Xi, defined accordingly. For anodic current 
Xu = Xcu + Xou + Ku, while for cathodic current, = Xcu + Xqu + Xau- Carrying out the various 
integrations leads to the result employed in the main article. 

i2 1 /' , NN /X Rewz 



la = e{nr) I Var T ^ / de{l - f{e))pa{e)-^ (A-21) 
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Figure 1: Comparison of density of states of the adsorbate for weakly coupled regime at low 
{9 = 0.1) and high coverage factor {6 = 0.9). he values of parameters (in eV) are as follows: 
E; = 0.6, E^^(O) = 0.2, = 0.4 and v = 0.5 eV. 
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Dos of adsorabate in strong coupling limit 
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Figure 2: Comparison of density of states of adsorbates for strong coupling regime at low 
and high coverage factor. The values of the various parameters employed (in eV) are as 
follows: El = 1.0, E^^(O) = 0.25, E^(0) = 0.75, An = 1.5, = 1.5, = 4.5, v = 2.0 
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Figure 3: Plots showing the density of states for redox, adsorbate and the Fermi distribution 
for anodic current under zero overpotential. The weakly coupled regime and low coverage of 9 
= 0.3 is considered here .The values of parameters (in eV) are as follows: = 0.6, E^^{0) = 
0.2, = 0.4 and v = 0.5 eV . 
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Figure 4: Plots showing the density of states for redox, adsorbate and Fermi distribution for 
cathodic current at zero overpotentiaL The values of parameters are same as in [3] 
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Figure 5: Plots showing the density of states for redox, adsorbate and the Fermi distribution 
for anodic current under zero overpotential. The strongly coupled regime and low coverage of 
6 = 0.3 is considered here .The values of parameters (in eV) are as follows: = 1.0, E^^{0) = 
0.25, E'^ = 0.75 and v = 2.0 eV . 
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Figure 6: Plots showing the density of states for redox, adsorbate and Fermi distribution for 
cathodic current at zero overpotential. The values of parameters are same as in [5] 
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Figure 7: Plots showing the variation of A(/) with respect to 9 the coverage factor. The values 
of re-organisation energies employed were same in both the curves. = 0.6 eV, Ea{0) = 
0.4 eV, EariO) = 0.2 eV 
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Figure 8: anodic current vs t] for a = 0.3. The values of the various parameters employed (in 
eV) are as follows: E; = 1.0,E^^(0) = 0.25, E^(0) = 0.75, A|| = 1.5, = 1.5, fi = 4.5, t; = 
2.0 



23 




Figure 9: anodic current vs rj for 6 = 0.1 in the weak coupled regime. The values of 
parameters (in eV) are as follows: E'!^ = 0.6, E^^{0) = 0.2, = 0.4 and v = 0.5 eV. 
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Figure 10: anodic current vs rj for 6 = 0.3 in weak coupled regime. The values of parameters 
(in eV) are as follows: E; = 0.6, ^^^(0) = 0.2, E'^ = 0.4 and v = 0.5 eV. 
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Figure 11: anodic current vs rj for 6 = 0.7. The values of parameters (in eV) are as follows: 
= 0.6, El^{0) = 0.2, = 0.4 and v = 0.5 eV. 
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Figure 12: anodic current vs 1] for 6 = 0.9. The values of parameters (in eV) are as follows: 
E; = 0.6, E^^(O) = 0.2, = 0.4 and v = 0.5 eV. 



27 




Figure 13: anodic current vs rj for 6 = 0.1. The values of the various parameters employed 
(in eV) are as follows: E; = 1.0,E^^(0) = 0.25, E^(0) = 0.75, A|| = 1.5, = 1.5, = 
4.5, t; = 2.0 
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Figure 14: anodic current vs rj for 9 = 0.3. The values of the various parameters employed ( 
eV) are as follows: E; = 1.0,E^^(0) = 0.25, E^(0) = 0.75, A|| = 1.5, = 1.5, = 4.5, t; 
2.0 
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Figure 15: anodic current vs rj for 6 = 0.7. The values of the various parameters employed 
(in eV) are as follows: E; = 1.0,^^^(0) = 0.25, E^(0) = 0.75, A|| = 1.5, = 1.5, = 
4.5, t; = 2.0 
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Figure 16: anodic current vs rj for 6 = 0.9. The values of the various parameters employed 
(in eV) are as follows: E; = 1.0,^^^(0) = 0.25, E^(0) = 0.75, A|| = 1.5, = 1.5, = 
4.5, t; = 2.0 
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Comparison of I ( r| = 0) for weakly and strongly coupled regimes 
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Figure 17: Plots showing the equilibrium current at zero overpotential /q vs 9 for strong and 
weak coupled regime. The values of re-organisation energies were selected be the same for 
both the curves, Er = 0.6 eV, ^^(0) = 0.4 eV, ^^^.(0) = 0.2 eV 
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Figure 18: Comparison of cathodic current obtained from corresponding expression for ca- 
thodic current and from empirical relation of anodic current, Ic = exp[—(3if]\Ia- The specific 
case selected is for 9 = 0.1 and a = 0.5 in the weakly coupled regime. The values of param- 
eters (in eV) are as follows: E; = 0.6, E^^(O) = 0.2, = 0.4 and v = 0.5 eV. 
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